On elliptical soft colloids in smectic-C films 
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We investigate theoretically the elliptical shapes of soft colloids in freely standing smectic C 
films, that have been reported recently. The colloids favour parallel alignment of the liquid crystal 
molecules at their surfaces and, for sufficiently strong anchoring, will generate a pair of defects at 
the poles of the colloidal particles. The elastic free energy of the liquid crystal matrix will, in turn, 
affect the shape of the colloids. In this study we will focus on elliptical soft colloids and determine 
how their equilibrium shapes depend on the elastic constants of the liquid crystal, the anchoring 
strength, the surface tension and the size of the colloids. A shape diagram is obtained analytically, 
by minimizing the Frank elastic free energy, in the limit of small eccentricities. The analytical results 
are verified, and generalized to arbitrary eccentricities, by numerical minimization of an appropriate 
Landau free energy. The latter is required for an adequate description of the topological defects 
when the liquid crystal correlation length is comparable to the size of the colloidal particles. 

PACS numbers: 61.30.Dk, 61.30.Gd, 61.72.Qq, 82.70.-y 
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I. INTRODUCTION 

Owing to their intriguing and complex behaviour, col- 
loidal dispersions in liquid crystals have been the sub- 
ject of numerous studies in recent years 1\. The be- 
haviour of these inverted emulsions depends upon (i) the 
elastic constants of the liquid crystal, (ii) the size and 
shape of the colloidal particles, (iii) the surface tension 
and (iv) the boundary conditions at the surface of the 
container. These contributions lead to highly anisotropic 
long-ranged colloidal interactions |3jQ1 ^^^^ result in a va- 
riety of self-organised colloidal structures, such as linear 
chains A,3.,G,T_,8i, periodic lattices 8, 9, 101 , anisotropic 
clusters [il|, and cellular structures [l^l stabilized, in 
general, by the presence of topological defects. 

The competition of elastic and surface ener gies may 
also determine the shape of nematic droplets \i3L Il4| . 
Indeed, for sufficiently large parallel anchoring, nematic 
droplets were found to exhibit sharp ends 15] . These 
shapes are known as factoids and the tactoidal shape is 
more pronounced for smaller droplets. As the size of the 
droplet increases, the contribution of the isotropic surface 
tension also increases and, as a result, the shape of the 
droplets becomes spheroidal .Igj. Recent studies on the 
shape of isotropic soft colloids in a nematic liquid crystal 
illustrate this behaviour and reveal a way of controlling 
the shape of the droplets using surfactants [l3| ■ 

Freely standing smectic films (Sm FSF's) are two- 
dimensional systems obtained by spreading a smectic liq- 
uid crystal through a hole drilled into a substrate. These 
films are useful systems to study surface interactions, ef- 
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fects of reduced dimensionality on liquid crystal phase 
transitions [igj £^rid interactions between soft colloids 
1^ Q, IM uM- ^oft colloids or inclusions are created by 
(i) rapid reduction of the film area, (ii) blowing across 
the film to pull material from the meniscus region or 
(iii) heating the system near the smectic-to-isotropic or 
smectic-to- nematic phase transition |2flj. 

The thickness of Sm FSF's ranges from a few thou- 
sand down to two smectic layers. Recent studies on 
Sm-C* FSFs at the Sm-C*-cholesteric (N*) phase transi- 
tion, revealed the existence of an intermediate thickness 
where N* soft colloids are nucleated ^^. For thin films 
{Af < 24; Af is the number of layers) the phase transi- 
tion is driven by layer-by-layer thinning, while for thick 
films {Af > 400) the N* order is stabilised and appears 
through the nucleation of A^* fingers. 

Experimental studies on the ordering of soft colloids in 
Sm-C* FSF report that the molecules in the liquid crys- 
tal matrix are aligned perpendicular to the colloidal sur- 
faces, nucleating one hyperbolic topological defect close 
to the colloidal surfaces [q|. By contrast, experiments on 
Sm-C FSF report that the liquid crystal molecules are 
aligned parallel to the colloidal surfaces T, 8, 10]. Under 
these conditions, a pair of defects may be nucleated at the 
poles of the colloidal particles which, in turn, affect the 
colloidal shape as reported by Cluzeau et al jTj]. Large 
colloids are found to be almost circular, while smaller 
ones are elliptical, with the long axis oriented along the 
Sm-C director, c. 

Here we carry out a theoretical investigation of the 
shape of soft colloids. We restrict our study to the ellip- 
tical shapes observed in the experiment referred to above 
13 and will not consider factoids. This article is orga- 
nized as follows. In section ^ we describe the Frank free 
energy for Sm-C FSF's. In Section ITTll we calculate an- 



alytically the configuration of tfie c-director for elliptical 
soft colloids, in the limit of small eccentricities. These 
results provide a complete (albeit approximate) descrip- 
tion of the dependence of the equilibrium shape (optimal 
eccentricity) on material parameters, that will be used 
to guide the subsequent numerical study. In Section llVI 
we use a Landau free energy to investigate numerically 
the same problem and compare the results obtained us- 
ing both approaches. Finally, in Section^we summarize 
our conclusions. 



II. FREE ENERGY 

The free energy of a Sm FSF takes into account the 
deformations of the in-layer molecular alignment, or ori- 
entational order, represented by a two dimensional vector 
field c = (cos 6*, sin 6*). It is given by 

T^ \ [ dS [lCs{V ■ cf + ICb{V X c)^] 

+ ^^rfz(7+Ysin2(r-0)), (1) 

where ICs and JCb are the splay and bend elastic con- 
stants, and 7 and W are the isotropic and anisotropic 
line tensions, respectively. 9'^ is the preferred orientation 
at the boundary dfl. In general, JCs and /Cf, are different 
but, for simplicity, we will use the one-elastic constant 
approximation, /Cs = /Cf, = /C. The bulk elastic free en- 
ergy density is then given by /C(V6')^/2. 

Minimization of Eq.lQJ with respect to 0{x,y) yields 



9 = in SI, (2) 

) = indn, (3) 

where v is the unit vector normal to the boundary dfl. 



v-VB- ^ sin 2(6*^ 
2/C 



III. ELLIPTIC DEFORMATIONS 

Let us consider a soft colloid in a smectic layer (Fig. 
(P)). The shape of the colloid is not fixed but depends, 
in particular, on the line tension. If the isotropic line 
tension is large, the colloid is circular, or slightly dis- 
torted. For simplicity, we start by considering elliptical 
deformations with small eccentricities. 

The problem of finding the equilibrium shape, or the 
optimal eccentricity, is simplified by using elliptic coor- 
dinates 



X = TZe cosh u cos ( 
y = TZe sinh u sin 4 



(4) 



where TZe — eTZ{l — e^) ^1^ , the eccentricity of the colloid 
is e = 1/ coshu'* and u^ is the minimal value of w, at the 
boundary of the colloid. TZ = y/A/n, where A is the area 
of the soft colloid. 



(a) 



(b) 





FIG. 1: (Color online) Soft colloids in a smectic-C layer with 
correlation length ^ — 0.17^. The vectors represent the c- 
director, while the colored regions correspond to strong vari- 
ations of the order parameter and signal the surface defects, 
(a) u — WTZ/IC = 1.0, a — 'yTZ/K, = 10 and optimal eccen- 
tricity e* = 0.270. (b) a; = 10, cr = 100 and e* = 0.127. 



In what follows we assume that the colloidal boundary 
favours a parallel orientation of the in-plane liquid crys- 
tal director. This will depend on the eccentricity of the 
colloid and is given by 



arctan 



yr 



sin 2(f> 



1 — cos2(j) 



(5) 



The solution of the Laplace equation© is straightfor- 
ward and has the general form 



6l(w,0) 



m=0 



exp {—mu} (him sin m(j} + Vm cos mq 



(6) 

for uniform liquid crystal alignment at infinity. Note that 
9{u, (/)) satisfies the following symmetries: (i) 9(u, —(p) = 
—9{u,(j)) and (ii) 9{u,'k — (j)) — —9{u,(j))^ implying that 
Vm = and restricting the nonzero coefficients to those 
with even indices, respectively. Redefining the coeffi- 
cients lA2m — > Z-/2m exp{2?nu*}, we write 

oo 
9{u,ip)— y Z^2m CXp{ — 2771 (w — U*)} sin 27Tl(/). (7) 

m— 1 

The analytical solution for elliptical soft colloids, for 
e <C 1, is obtained using perturbation theory. We expand 



/(O) 



the coefficients U2m — ^2m "^ 
e^92 + e^9i 



c^U?^ 



■e4W^l^ + ..., ofthe 



(4) 
2m "^ ■=• ^'2m 

solution 6' = 6'o + e'^6'2 + e'*6'4 -I- ..., and generalize the 
method used for circular colloids in [23 ■ Assuming that 



U 



2m = —Po'/m, for circular soft colloids, 6*0, is given by 
Po exp {— 2(u — u*)} sin 2(^ 



arctan 



1 -poexp{-2(i 



')}cos2(/) 



(8) 



where po is obtained by minimizing the free energy 
Eq.(^. This yields a second order algebraic equation 



pI 



+ (4/c^)po -1 = 0, 



(9) 



with a positive root po = {'^/^) \/l + (^/2)^ — 1 
where lo = WR-/K. is a dimensionless ratio of the an- 
choring and elastic strengths. 



Suppose we can write the free energy density as f[9] = 
/o[^]+e^/2[0] + e'*/4[^] + ... A functional Taylor expansion 
of the free energy about the solution for a circular colloid, 
6*0, yields 



J^[9] = J d^^ifoieo] 



e'hm+e^hm) 



=2 <5/2 



59(x) 



i/d2x'/d2xA0(x) 



s'h 



se(x)se{x') 



Oo 



A0(x') + ..., (10) 



where A6 — 6 — Oq = e 02 + e 0^ + ..., and the expansion 
is valid to the fourth order in the eccentricity. This is 
simplified using the equilibrium condition for the circular 
case, (5/o/(50|g = 0, which implies that 02 and 0^ are 
not present in the terms of order e^ and e^, respectively, 
yielding for the free energy 



T[0]^Tom 



e^T2\ 



e^TA 



(11) 



The transition from circular to elliptical shape is con- 
trolled by the sign of ^2, which depends only on ^o- If 
this term is positive, the soft colloid is circular; otherwise 
it is elliptical. To determine the optimal eccentricity, we 
begin by minimizing the free energy fEa. ljllll ) with re- 
spect to 02- The solution satisfies the linear equation 
5!Fi/502 — 0, and may be written as 



'«-/ 



d' 



S'fo 



60{x)S0{:ic') 



Sh 



(x') 



(12) 



eo 



The general solution is a linear combination of the 
solutions of the Laplace equation, in terms of which 
the linear equation becomes a matrix equation for the 
set of coefficients Uj^^^ ^ind Eq. (|12|l an (infinite) ma- 
trix inversion. The optimal eccentricity is then, e* ~ 

V-^2 [001/2.^4 [00, 02]. 

In order to determine To, T2 and J^4 in Eq. (fTT|l . we 
proceed by noting that the bulk elastic free energy may 
be integrated by parts, yielding 



^k//C 



1 r' 

2 Jo 
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dcj)0 — 

ou 


on 
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Jo 


',02 90-, 
^^^2 du 



-i^ > (13) 

where — /Cfobo] — — T^/Clog (1 — Po) is the elastic free 
energy of a circular soft colloid. Terms of order e^ and 
e^, that depend on 02 and 04 respectively, were neglected 
since they do not contribute to the total free energy. 

The term corresponding to the isotropic line tension 
7 is independent of 0. Its expansion in the eccentricity 
yields 



T^jK, 



-i/JC / dl 



2vra ( 1 + -e 



(14) 



where we used dl/dcf) — TZe\/l — e^ cos^ 0/e. a = •jTZ/JC 
is the ratio of the isotropic line tension and the elastic 
strength. 



The expansion of the contribution of the anisotropic 
line tension or anchoring strength, W, is obtained by 
expanding the preferred orientation (Eq.©), 0'^ — 0q + 
e^0| -I- ... After some tedious algebraic and trigonometric 



manipulations, we obtain 


T /r ^/^ f 
Tw/JC = / ° 

^ Jan 


^ 2 


So\po] 4 



dl sin' 



'S2[po] 



20?, 



Si[po] 



,2{0l 



2-rr 



i^ COS 20 sin 2 (0g 



(15) 



where /Ca;5o[po]/2 = 7r/Ccj(l —po)/2 is the anchoring en- 
ergy for a circular soft colloid, 52[po] = — 37r(l — Pq)/8 
and 54[po] = 7r(-18 -h 5po + 24pg - 3pi^)/128. 

By adding Eqs.JEI), (|I3I> and (jTHl) we obtain the free 
energy, T[0] ~ Ta[pa] + e'^T2[po] + e**J^4[po, 62], where 



To[po]/K. 
^2[po\/K, 



-£o[p0]+27rcr 

ujS2[po\/2 



-ujSo[po\/2, 



(16) 
(17) 



and 



Ti[pa,02]/K. = -Silpo] 



— TTcr 

32 2 



502 

du 



on 



r2Tr 



2(02 -20^) cos 2(0^5-00) 

UJ f^ 

+ - d(P02Cos2(l>sm2{0'o-0o). (18) 

8 Jo 

Inspection of ^-"2 [po] reveals that the shape transition is 
controlled by the dimensionless quantity, lo. It is clear, 
from the definition of po, that ^2 is negative for any finite 
value of LO and thus the equilibrium shape is elliptical, 
except when to vanishes (weak anchoring). 

In order to calculate 02 , we approximate the general so- 
lution by a finite number of terms in ((TJ, corresponding 

(2) 
to a small number of coefficients Zi2,„ . Explicit calcula- 

(21 (2) 

tions reveal that the first two coefficients, U2 and Ul , 
are sufficient to yield 02 and the eccentricity e* with an 
accuracy better than 1% and 0.6%, respectively. 

Figure Q illustrates the dependence of U^ and U^ 
on the reduced anchoring strength u>. For weak anchor- 
ing, as UJ vanishes, we find 



^2 - ^'^^ 



7/(2) _ (^ 

"4 - ^64 



(w<l). 



(19) 



while for strong anchoring, 

/(2) _ 1/1 1\ 1 
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2aj2 



u^ 
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{If 



{U > 1). 



(20) 
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FIG. 2: U\ and U\ as functions of the reduced anchoring 
strength u = WIZ/K.. Higher order coefficients lead to small 
corrections (less than 0.6%) to the optimal eccentricity e*. 



FIG. 4: Shape diagram: lines of constant eccentricity in the 
uj-a plane. The circle and the ellipse schematically shown 
occur for e = and e = 0.9, respectively. 
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FIG. 3: Optimal eccentricity e* as a function of a = "jTZ/K, 
at fixed values of c^ = WTl/K. = 0.1, 1.0, 10, oo. 



In Fig. © we plot the optimal eccentricity e* as a 
function of a^ at fixed values oi oj — 0.1, 1.0, 10, oo. As 
expected, the stabiUty of the elliptical shape decreases as 
(J increases. We note, however, that the rate at which 
the eccentricity e* decreases as a increases depends on 
Lo. In fact, as oj increases a larger isotropic tension, a, 
is required to reduce the eccentricity. In the limit of 
strong anchoring, the liquid crystal molecules are parallel 
to the surface of the soft colloid, but at the same time, 
are parallel to each other, to minimize the elastic free 
energy. These conditions imply a strong eccentricity, or 
a pronounced elliptical shape, of the soft colloid. The 
isotropic line tension, on the other hand, favours circular 
soft colloids, minimizing the perimeter and reducing the 
eccentricity. 

For weak anchoring, the optimal eccentricity has the 
asymptotic behaviour 



This expansion clearly reveals a bifurcation at w = 0. At 
zero anchoring, the liquid crystal is not aligned by the 
colloidal boundary, and the shape of the soft colloid is 
circular. For positive values of w, the soft colloid becomes 
elliptical. For strong anchoring, the optimal eccentricity 
,e*, has the asymptotic behaviour 



2 73(19 + 12cr) 1 

v/a-ll/3~ 2(3^-11)3/2^ 



{uj > 1). (22) 



(w/a)i/2 + 3/8(c^/a)^/2 (^ « i). 



(21) 



When ijj ^ oo the colloidal shape is strongly affected 
by the elastic deformation of the liquid crystal matrix 
and the eccentricity is determined by the competition 
between the isotropic line tension and the elastic free en- 
ergy. When a is large, the eccentricity is small, and the 
soft colloid is almost circular. As a decreases, the soft 
colloid becomes more elliptical and, eventually, this per- 
turbation theory breaks down. We note, however, that 
the absence of solutions of (|22|1 for a < 11/3 may indicate 
the stability of other colloidal shapes, such as the factoid, 
an elongated shape with two singularities or cusps at its 
extremities. 

Figure Q) illustrates the shape "diagram" of the soft 
colloid depicting lines of constant eccentricity in the Lo-a 
plane. These curves exhibit two regimes. For w < 1, the 
slope of the constant eccentricity lines indicates that a is 
proportional to uj. For larger values of w, however, the 
lines are nearly flat, indicating a dependence on a only. 
The lower line on the diagram corresponds to a "critical" 
value of (Tc where the eccentricity of the colloid is e* = 1 
and the elliptical shape becomes unstable. 

The shape diagram reveals that colloids with approx- 
imately circular shapes (e — *■ 0) occur when a is 2-3 or- 
ders of magnitude larger than uj, in the weak anchoring 
regime. In the strong anchoring regime, nearly circular 
colloids occur at large values of a with a dependence on 
UJ that is less marked. This means that the shape of the 
colloids depends not only on the ratio of the isotropic 



to anisotropic line tensions o Jlo^ but also on the director 
configuration, as expected. 



IV. NUMERICAL SOLUTIONS 

The Frank elastic free energy used in the previous sec- 
tion (Eq. (^) describes the general features of soft col- 
loids in liquid crystals. However, it does not describe 
adequately the topological defects that may occur in the 
strong anchoring regime. The structure of these defects is 
relevant for colloids with sizes comparable to the liquid 
crystal correlation length. Furthermore, the analytical 
solutions described above are based on a perturbation 
theory that is valid for small eccentricities only. 

For completeness, we consider a Landau-like free en- 
ergy, where the orientational order parameter of the Sm- 
C phase, ^ = |'0| exp [iff), includes both the tilt angle |-i/;| 
and the azimuthal orientation Q. The Landau free energy 
may be understood as a Taylor expansion in terms of the 
invariants of ?/> and of its deformations or spatial deriva- 
tives dii\). For an isotropic Sm-C film (with no preferred 
orientation within the layers), there is a single quadratic 
and a single quartic invariant of V'- In what concerns the 
derivatives of ?/>, however, several quadratic invariants 
exist but they are reduced to a single term proportional 
to |V'!/'p = l^j^V'P + |f^yV'p5 in the one-elastic constant 
approximation. 

We use this notation to rewrite the surface free energy, 
Eq.(^, denoting by V'^ the preferred order parameter at 
the boundary of the colloid. The surface free energy writ- 
ten in terms of an order parameter such as ?A has to ac- 
count for (i) the coupling between the orientation of the 
smectic molecules and the surface and (m) the coupling 
between the tilt angle j-^l and the one imposed by the 
surface |-0*|, resulting in a form that is slightly different 
form that of Eq. (QJ . The total free energy is 
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dS 
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M 
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C. 



dl 



dn 
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1 



jV,|4 + _|vv;| 



(23) 



where Re{} denotes the real part and ip* is the com- 
plex conjugate of tJj. Using ipbuik — \/A/B, for the order 
parameter of the uniform Sm-C, and making the trans- 
formations V' = V'/(V2|'06«;fc|)i ^ = x/TZ, we obtain the 
reduced free energy T 



Ti{An^\'ilJbuik?) 



T^ dS 
hi 



|V'P(|V'P-1) + £|V7/^P 



dl 



dil 



I -Re 






(24) 



where e = [£^/Rf and ^ = \/C/A = yjK./AI\i^buik\ is 
the correlation length of the Sm-C, the typical length of 
a defect. The Frank free energy (Eq. (Q), is recovered 




FIG. 5: Reduced free energy T as a. function of the eccentric- 
ity for e = 10~*, uJ = 10 and a = 6. 



when e -^ 0. For simplicity, we have considered \ip^\ = 

\'4^buik\- 

The size of the colloids may be 10 to 100 times the 
size of the defects, i.e., e <C 1. The major difficulty in 
the numerical solution stems exactly from these different 
length scales. We use finite elements with adaptive mesh- 
ing, as described in Ref . [l9l| . to minimize Eg. 1)24(1 . A first 
triangulation respecting the predefined boundaries is con- 
structed. The order parameter ^ is given at the vertices 
of the mesh and linearly interpolated within each trian- 
gle. The free energy is then minimized using standard 
methods. The variation of the solution at each iteration 
is used to generate a new mesh. In typical calculations 
convergence is obtained after two mesh adaptations, cor- 
responding to final meshes with 10'* points, spanning a 



region of 307?. x 307?., and minimal mesh sizes of 10^^7?., 
close to the defects. The free energy is obtained with a 
relative accuracy of 10"'^%. 

The optimal eccentricity is obtained as follows. We 
start by fixing the set of material parameters {e, to, a}. 
For a given value of the eccentricity e the free energy 
of Ea. (|24|) is minimized. The eccentricity is incremented 
by 5e and the minimization is repeated. The free en- 
ergy profile, as a function of e, exhibits a well defined 
minimum (in most cases) , allowing the calculation of the 
optimal eccentricity e*(Fig. |(SJ|). For the smallest ec- 
centricities e*, however, the variation of the free energy 
near the minimum !F{e*) is of the order of the relative 
numerical accuracy. This renders the method less accu- 
rate in the limit of circular colloids (see Fig. ^). After 
determining e*, the parameter a is incremented and the 
free energy profile is calculated. 

In the previous section we calculated the optimal ec- 
centricity, e*, analytically. Since the results are based 
on perturbation theory, the numerical results are not ex- 
pected to match the analytical ones in the limit e — > 0, 
for any value of e* . Figure 10 illustrates how the optimal 
eccentricity depends on e for a) w = 1 and b) lo = 10, and 
a — 10. Both curves approach a constant value of e*, as 




100 



FIG. 6: Optimal eccentricities as a function of e for a) w = 1 
and b) tj = 10, and a — 10 
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FIG. 7: Optimal eccentricity e* as a function of a at constant 
Lu — 1, 10, 100 and e = 10~^, 10~*. Lines and symbols 
represent the analytical and numerical results, respectively. 



e —^ and the numerical results approach the analytical 
values, when e — > 0, for small eccentricities. 

In Fig. lO we plot the numerical results together with 
the analytical curves for the optimal eccentricity as a 
function of a, at constant w = 1, 10, 100 and e = 10~^, 
10^"*. The numerical results confirm that the optimal ec- 
centricity is a strictly decreasing function of the line ten- 
sion a. Moreover, the decrease in the slope of e* versus 
a is more pronounced in the presence of surface defects 
{uj > 10). This slope does not change significantly with 
s (at constant uj). 

The numerical results for lo — 100 cross the corre- 
sponding analytical curve, which in view of the results 
for lower w, may seem surprising. This crossover was 
observed for anchoring strengths, in the range 10 < a; < 
100, and is due to the fact that the critical value CTc, where 
the elliptical shape becomes unstable, (e* -^ 1), does not 
behave as predicted by the analytical theory. In fact, 
while the analytical approach predicts that ctc/uj varies 
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FIG. 8; Shape diagram. Lines and symbols denote equilib- 
rium shapes at constant eccentricity, in the uu — a plane. The 
lines were obtained analytically while the symbols are the re- 
sult of a numerical minimization, for systems with e — 10~^. 



considerably, the numerical results suggest a nearly con- 
stant value, ac/u! = 0.505 ± 0.005. It is likely, however, 
that this limit will not be observed in experiments since 
for most materials a/ui ^ 1 |24l |. 

In Fig. JHJ we plot the numerical and the analyti- 
cal shape diagrams for several values of the eccentricity 
(e = 0.1,0.2, 0.3, 0.4, 0.5). The numerical shape diagram 
was obtained for e = 10^^. Results for e = 10""' yield 
similar curves, that are shifted by a small positive value 
of a. The results of the numerical calculations reveal the 
existence of two distinct regimes. For weak anchoring, 
UJ < 10, circular colloids may occur if a is 2-3 orders 
of magnitude larger than to, in line with the analytical 
results. For strong anchoring, lu > 10, the presence of 
the defects leads to a behavior that differs from the an- 
alytical predictions. In fact, the latter suggests that a is 
nearly independent of w, while the numerical results in- 
dicate that a increases as lo increases. In the limit e = 1, 
the results of the numerical calculations suggest a linear 
behavior, ar ~ 0.505a;. 



V. CONCLUSIONS 

We have studied the shape of soft colloids in Sm-C 
FSF's. We assume that the colloids may adopt one of two 
shapes (i) elliptical with eccentricity e and (ii) circular 
(e = 0) and using the Frank elastic free energy obtained 
an approximate solution of the c-director field valid for 
small eccentricities. This provides insight on how the 
eccentricity depends on the material parameters, such as 
the dimensionless quantities lo and a. The analysis also 
shows that soft colloids adopt non-circular shapes for any 
finite anchoring strength, lo. We have found a critical 
value of (Tc where the eccentricity is e* = 1, indicating 
that the elliptical shape is no longer stable. 



The analytical results reveal that the size of the col- 
loid TZ plays a very important role. For weak anchoring 
conditions, WTZ < /C, colloids with approximately circu- 
lar shapes may occur if 7 is 2-3 orders of magnitude 
larger than W. For strong anchoring, circular shapes 
occur when TZ > /C/7. In all cases, elliptical shapes are 
favoured for sufficiently small colloids. 

Finally, we compared the analytical results with the 
results of a numerical minimization of a Landau-like free 
energy. We have proposed an expression for the surface 
free energy that takes into account the coupling between 
the orientational order parameter of the Sm-C phase and 
the order favored at the boundary of the colloid. The 
numerical results validate and extend the results of the 
analytical calculations. The main difference concerns the 



curves of constant eccentricity depicted in the shape di- 
agram of Fig.JHJ, in the strong anchoring regime. While 
the analytical approach predicts curves that are almost 
constant for large uj, the numerical results indicate that 
a continues to increase as to increases and in the limit 
e = 1, ac, is found to vary linearly with lo. 
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